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Abstract 

We study the energetics of the accretion-induced outflow and then plau- 
sible jet around black holes/compact objects using a newly developed disc- 
outflow coupled model. Inter-connecting dynamics of outflow and accretion 
essentially upholds the conservation laws. The energetics depend strongly 
on the viscosity parameter a and the cooling factor / which exhibit several 
interesting features. The bolometric luminosities of ultra-luminous X-ray 
binaries (e.g. SS433) and family of highly luminous AGNs and quasars can 
be reproduced by the model under the super-Eddington accretion flows. Un- 
der appropriate conditions, low-luminous AGNs (e.g. Sagittarius A*) also 
fit reasonably well with the luminosity corresponding to a sub-Eddington 
accretion flow with / — > 1. 
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1 Introduction 



Extremely high resolution observations of the powerful extragalactic dou- 
ble radio sources HBlandford & Rees 19741 [Begelman, Blandford & Rees 1984 



Ferrari 19981 IMirabel 200311 reveal that they are formed by well-coUimated out- 
flows or jets that continuously emerge from the nuclear region of the host active 
galaxies (AGNs) or quasars, believed to harbor supermassive black holes. Simi- 
larly micro-quasars p Mirabel & Rodriguez 1994] Mirabel & Rodriguez 1998 1 dis- 
covered in recent times reveal that outflows are generated from stellar mass black 
holes (or black holes candidates). The well-coUimated outflow in SS433 ob- 
served for several decades, which is our galactic, persistent, super-critical ac- 
cretor, is an well known evidence of cosmic jet JMargon 1984 1. Further, highly 
coUimated line jets are seen in young stellar objects HMundt 198511 . Relativis- 
tic jets are also observed in neutron star low mass X-ray binaries (LMXBs) 
flMigliari & Fender 2006p . 

The outflows/jets extract matter, energy and momentum from the accretion 
disc that forms around the compact object, and thus the dynamics of outflow lead- 
ing to jet is intrinsically coupled with the accretion dynamics through the conser- 
vation laws. Also several observations (precisely the simultaneous observations of 
disc and jet; see e.g. Ghosh & Mukhopadhyay 2009, and references therein) reveal 
that accretion processes and outflows are strongly correlated and they eventually 
control the accretion process, precisely, the accretion dynamics in the vicinity of 
the central star. The relativistic outflowing matter, in the case of quasars or micro- 
quasars, should indeed come only from the inner region of the accretion disc. 
This is particularly suggestive as the quasars or the micro-quasars do not have 
an atmosphere of their own. However, most of the models of accretion disc and 
outflow/jet have been evolved separately, considering these two to be apparently 
dissimilar objects. The outflow/jet models have been evolved over the years, from 
speculative ideas such as de Laval nozzles IBlandford & Rees 19741 to electrody- 
namic acceleration model piandford & Znajek 1977[ , centrifugally driven out- 
flows [ [Blandford & Payne 1982llPudritz & Norman 19861 . etc. Moreover, MHD 
simulations of outflow/jet have been performed both in non-relativistic as well as 
relativistic limits, mostly in the Keplerian paradigm (e.g. Hawley & Balbus 2002, 
Mizuno et al. 2006, Hawley & Krolik 2006, and references therein) to see how 
the matter gets deflected from the equatorial plane. Nevertheless, the definitive 
understanding of the origin of outflows/jets is sill unknown. 

The radiatively driven outflow or jet whose origin is better understood, can 
be envisaged when the accretion disc is highly radiation pressure dominated or 
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precisely "radiation trapped". This is likely to occur when the accretion rate 
is super-Eddington or super-critical HLovelace et al. 19941 [Begelman et al. 20061 
IFabrika 20041 [Ghosh & Mukhopadhyay 2009 1 as in ultra luminous X-ray (ULX) 
sources such as SS433 (with luminosity ~ 10"^" erg/s or so) [IFabrika 200411 : a 
prototype of ULXs in external galaxies having relativistic jets. The super-critical 
flows are optically thick and have a strong advective component ^Lipunova 1999 



Ohsuga et al. 2005| . The unusual signatures of ULXs, compared to their counter- 
parts, lie in their extreme high luminosity (L^o/ ~ iQ^^-'^^ erg/s), strong spectral 
variability and their association with the actively star forming regions, that make 
their nature a subject of controversy IFa brika 20041 Mushotzky 2004 1. Recently, 
an ultra-luminous accretion disc with a high kinetic luminosity radio jet has been 
discovered in quasar PKS 0743-67 fPunsly & Tingay 2005 1, which indicates that 
the ultra-luminous sources may be extended to quasars or even AGNs. On the 
other extreme end, the under-luminous AGNs and quasars (e.g. Sagittarius A*) had 
been described by the advection dominated accretion flow model (in short ADAF 
[ |Narayan & Yi 1994| ), where the flow is substantially sub-critical/sub-Eddington. 
The model hypothesized a plausible emergence of strong outflows/jets, which lead 
to reveal a strong inter-connection between the outflow and advection. 

Very few models exist which simultaneously deal with the accretion and out- 
flow/jet dynamics on the same platform, both in analytical regime as well as 
through numerical simulations (see Ghosh & Mukhopadhyay 2009, and refer- 
ences therein). The difficulty of simultaneous simulation of the disc and outflow 
arises due to the diff"erence in time scales between accretion and outflow. Also 
the results strongly depend on the initial conditions [Ustyugova et al. 19991 • The 
2.5-dimensional disc-outflow coupled model given by Ghosh & Mukhopadhyay 
(2009) has been formulated by incorporating explicit information of the outflow 
in a fully analytical regime through a self-similar approach in a general advective 
paradigm, by upholding the conservation laws. It has been further shown that 
with mass, energy and momentum conservation, and a few scaling arguments, the 
correlation of the jet with the disc could be successfully modeled. The authors 
solved a complete set of partial diff"erential fluid equations (with the variation of 
flow parameters in both radial and vertical direction) without assuming the vertical 
hydrostatic equilibrium, with the explicit inclusion of the vertical velocity (repre- 
senting outflow) and all the relevant components of the stress tensor apart from the 
usual W,-^, from the first principle. They used their solutions to study two extreme 
cases of the geometrically thick advective accretion flows: super-critical and high 
sub-critical, which are more probable regimes of the strong outflows and jets. 
They found that the flow parameters of the accretion-induced outflow/jet strongly 
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depend on Shakura & Sunyaev viscosity parameter a JShakura & Sunyaev 1973p 
and cooling factor /. 

In the present paper, we use and extend the above mentioned work 



[Ghosh & Mukhopadhyay 2009 1 to study the detailed energetics of the accretion- 
induced outflow and then plausible jet. In order to do that, we stick to two ex- 
treme regimes of the black hole accretion: super-critical and high sub-critical, 
to describe ultra/highly luminous and under-luminous sources respectively. For 
obvious reasons, we do not repeat the calculation and the formulation of the disc- 
outflow model given earlier [Ghosh & Mukhopadhyay 2009 1, but will recall them 
appropriately. 

We arrange the paper in the following manner. In the next section, we formu- 
late the equations for the energetics of the accretion-induced outflow. In §3, we 
study the properties of the disc-outflow energetics. Finally we end in §4 with a 
discussion and implications. 



2 Energetics of the accretion-induced outflow 

In order to study the energetics of the outflow, we need to compute the mass 
outflow rate. Blandford & Begelman (1999) generalized the ADAF solu- 



tion [Narayan & Yi 1994 1 by including wind and outflow, assuming the mass 
inflow rate to be proportional to r^, where < /> < 1 (r be the ra- 
dial coordinate). Chakrabarti and his collaborators P Das & Chakrabarti 1999[ 
IChakrabarti 199911 also attempted to calculate the mass outflow rate from a 
disc involving a shock without including the vertical flow explicitly. In the 
present paper, we derive this self-consistently. Integrating the continuity equation 



[Ghosh & Mukhopadhyay 2009 1 vertically about the equatorial plane we obtain 



£ r 

dr j/,„ 



Anrpvrdz + Anr{p{h)v,{h) - p{ho)v.{ho)) = 0, (1) 

ho 



where and are the radial and the vertical velocity respectively, h{r) is an 
arbitrary disc scale height and /iq the minimum (finite) height of the disc-outflow 
system, which is diff"erent for super-Eddington and sub-Eddington accretion flows, 
z be the vertical coordinate. Here we deliberately focus on the disc-coronal region 
(sub-Keplerian halo), from where the mass loss takes place in the form of wind 
or outflow, and hence the solution at equatorial plane looses significance. The 
integral in the first term of the left-hand side of eqn. ^ represents the disc mass 
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accretion rate given by 



Ma(r) = - f Anrpvrdz. (2) 

The second term then can be attributed to the rate of change of outflow (and then 
plausible jet) mass, given by 

dMj{r) dMair) 

— ^ = -Anr{p{h)vM - p(ho)v,(ho)) = f^. (3) 

dr ar 

Here, Mj{r) is analogous to the mass outflow rate represented as 



M 



4j(r) = - J 4nr(p(h)v,(h)-p(ho)v,(ho))dr + cj, (4) 

where the constant cj is to be determined by an appropriate boundary condition. 
Thus from eqn. ([3]) the total mass accretion rate can be written as 

M = Ma{r) + Mj{r). (5) 

Eqn. ^ entails that under the stationary condition the radial mass flux in the disc 
decreases as the inflowing matter approaches the central object, at the same rate at 
which the vertical mass flux increases to maintain a constant M which is exactly 
the net mass accretion rate at infinity. 

To compute the power of the outflow and then plausible jet, following previous 
work [ Ghosh & Mukhopadhyay 2009 J , we integrate the disc-energy conservation 
equation for an accretion-induced outflow over the disc scale height which yields 



d_ r 

dr Jha 



h 

AnrTrdz + Anr{TM - T.ih)) = 0. (6) 

ho 



Here, is the vertical component of the total energy flux given by 

= (y + + 0g)pv, - v^W^, + F„ (7) 



where = vj + + v^, (pc is the gravitational potential, the radiative flux 
from the disc surface, W^j. the ^z"' component of the stress tensor, y the ratio of 
the specific heats of the gas-radiation mixture. We thus obtain the power of the 
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outflow from eqn. which is the total power removed from the disc by the 
outflow, as 



Pj{r) = J AnriTM - T,{K))dr. (8) 
We then calculate the disc luminosity in presence of outflow (and jet) using disc- 



energy conservation equation described earlier [Ghosh & Mukhopadhyay 2009 1 



L = {l-f) J I QHnrd^ dr, (9) 

where is the total viscous heat generated in the disc. We now introduce 
three dimensionless parameters which correlate the disc and the outflow lead- 
ing to jet: qjm{r) = Mj(r)/M, qjpir) = Pj(r)/Mc'^ and qi = LjMc^. We do not 
describe here the generalized Bernoulli's number {Be), explicitly shown earlier 



[Ghosh & Mukhopadhyay 2009 1, which represents the total energy of the system. 



3 Properties of the energetics of disc and outflow 

Accretion by a black hole (or any other central star) is the primary source of 
the mass outflow and jet formed in the inner hot region of the disc. At the first 
instant we neglect the contribution of magnetic field which is likely to be the ori- 
gin of coUimation and acceleration of jet. However, doubts can be raised about 
the significance of the magnetic field for production, coUimation and accelera- 
tion of jets in ULXs and highly luminous AGNs and quasars, which are super- 
critical and highly radiation trapped systems (see Fabrika 2004, for details). We 
argue that the radiation pressure is likely to be the plausible reason for strong out- 
flows and then plausible jets in high mass accretion flows BLovelace et al. 1994|| . 
The inflowing matter is expected to be ejected and accelerated through a funnel 
like region PFukue 1987[ [Fabrika 200411 by the strong radiation pressure leading 
to strong outflow and then plausible jet from the disc. In early, based upon the 
geometrically thick accretion disc model, it was shown that very narrow and deep 
funnels are formed around the rotation axis of the accretion disc and most of 
the disc energy flux is radiated from the surface of the funnels. The radiation 
pressure may accelerate the outflowing mass in funnels up to relativistic veloci- 
ties in the form of two jets pointing in opposite directions as observed in SS433 
Paroszyriski, Abramowicz & Paczyhski 1980] [Abramowicz & Piran 19801 . In 
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the inner accretion disc, a few Schwarzschild radii away from the black hole, 
as the infall time scale is much smaller than the viscous time scale to transport 
angular momentum outside, matter has no time to loose angular momentum and 
attains a near constant value. Around this region the centrifugal force becomes 
comparable to the gravitational force, and thus the incoming matter slows down at 
this centrifugally dominated region and gets puffed up HChakrabarti 199911 . There- 
fore, a funnel type geometry is likely to form around the region and the radiation 
pressure blows up the matter through them. 

However, the magnetic field is likely to play the role in generating turbulence 



in the hot disc systems [Balbus & Hawley 1991 1. In fact, for critical or sub-critical 



accretion flows, the jet is likely to form due to the magnetic activity in the disc 
JBlandford & Payne 1982[|Camenzind 1986II. We do not intend to discuss here the 



possible reason of the emanation of the jet. Whatever the explanation/mechanism 
be, the potential dynamics of the accretion-induced outflow and then plausible jet 
always upholds the conservation laws which shows several interesting features. 
The dimensionless parameters qjmir), qjpir), qi and Be carry the information of 
the disc-outflow energetics. In the next two subsections we analyse their features 
for super-critical and sub-critical accretion flows respectively. 

3.1 Super-critical accretion 

To study the high mass accretion flows, we consider M ~ 3 x 10~'*MQ/yr corre- 
sponding to the central mass M ~ IOMq, which is a probable estimate of the ac- 
cretion rate of our own galactic super-critical accretor SS433 UFabrika 2004.1 . We 
express M in terms of critical Eddington unit so that, for our case, M ~ 10"* M^r- 
This makes the energetic profiles independent of mass of the central star and one 
can extend our model to the case of high or ultra-luminous AGNs and quasars. 
Figures la and lb depict the variation of mass outflow rate and power of the out- 
flow and then plausible jet respectively as functions of r. The nature of profiles 
is similar to that of vertical component of velocity v,, explained in detail earlier 



[Ghosh & Mukhopadhyay 2009 1. At / = 0.4, a substantial amount of matter and 
radiation is extracted in the form of strong outflow leading to jet and thus the 
power of the outflow is significantly high, which falls off rapidly for / > 0.4. 

The luminosity profile (Fig. Ic) exhibits an interesting feature. At small 
r, it increases slowly (although not very clear from the figure that is given in 
the logarithmic scale) with r and eventually attains a near constant value. At 
high a, the luminosity corresponds to the Eddington or the super-Eddington ac- 
cretion. However the system becomes sub-Eddington at low a ~ 0.01, even 
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if M ~ lO'^Mcr and / low. This extreme sensitivity is due to the fact that a 
directly controls the viscous heat generation. At / = 0.4, for a ~ 0.3, L 
reaches ~ lO'^^erg/s, resembling highly luminous and ultra-luminous XRBs (e.g. 
CXOU J095550.2+694047 UKaaret et al. 200 IL GRS 1915+105, GRO J1655-40 



[Mirabel & Rodriguez 1998 1). At / > 0.4, L decreases in orders of magnitude 
but still lies within the range of the high or ultra-luminous regime (as in the case 
of SS433). For a black hole of mass M ~ 10^ - WMq, it can be shown that L 
attains a value ~ 10"^^ - lO'^^ erg/s for a = 0.3 and / = 0.4. Such luminosity is 
commonly observed in highly luminous AGNs and ultra-luminous quasars (e.g. 
PKS 074367), possibly in ULIRs PGenzel et al. 199811 and narrow-line Seyfert 1 
galaxies (e.g. JMineshige et al. 2000| ). The nature of the variation of Be has been 



shown in Fig. ld,e, already explained by Ghosh & Mukhopadhyay (2009) with 
three-dimensional plots. 

The truncation of the curves attributes to the truncation of the disc due to evap- 
oration to corona. With the decrease of /, which corresponds to more increase of 
radiation pressure, the disc gets more radiation pressure dominated and the matter 
gets blown up due to the strong radiation pressure further out. At this situation, the 
disc tends to become more Keplerian and thus centrifugally dominated. Therefore 
the centrifugal dominated boundary layer is likely to form at a larger radius. As 
a result the outflow rate and then the outflow/jet power increase significantly and 
the disc gets truncated at a larger radius. 

In Fig. 2, we show a comparison of v^, qj^, qjp and qi between two values 
of 7 at the same set of a and /. With a small increase in the value of y from 
1 .4 to 1 .444, the above parameters decrease by an enormous margin. The abrupt 
change of the corresponding parameters reflects that the system is very sensitive to 
7 and then /3 which is the ratio of the gas pressure to total pressure Q. An increase 
of y from 1.4 to 1.444 corresponds to an increase of /3 from 0.3 to 0.5, which 
further corresponds to the conversion from radiation dominated to marginally gas 
dominated system. This signifies that the gas pressure in the system increases, 
which means that the radiation pressure of the system suddenly decreases. The 
decrease in the radiation pressure inhibits strong outflow [^. Hence the power of 
the outflow and mass outflow rate fall in order of magnitudes. The same reasons 
are applicable in order to explain the variation of luminosity profiles. 



'The relation between y and /3 is given in Ghosh & Mukhopadhyay 2009. 
^In super-critical, radiation pressure dominated regime, strong radiation pressure blows up the 
matter to form outflow and consequently plausible jet. 
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3.2 Sub-critical accretion 

Here we consider that the flow has a mass accretion rate M ~ 10'^ M^r correspond- 
ing to M ~ 10^ - IO^Mq, which renders the flow to be thick and gas pressure 
dominated. Figure 3 shows that the mass outflux and the power extracted from 
the disc by the outflow are much less compared to that of the high mass accretion 
flows. Even for a near extreme / ~ 0.7, 5^ is ~ 4 times larger for super-critical 
flows compared to the ideal case of an advection dominated (/ — > 1) sub-critical 
accretion flow. This signifies that outflows are more probable for super-Eddington 
accretion flows compared to that of sub-Eddington flows. For low a and / = 0.9, 
we obtain a luminosity L ~ lO^^ergl s corresponding to M ~ 10^ (Fig. 3c). 
For high a ~ 0.3, L increases by three orders of magnitude (not shown in the fig- 
ure) at the same /. Such low luminosities are observed in many under-luminous 
AGNs (e.g. Sagittarius A*; [IMahadevan 1998II ). 

4 Discussion 

Based on a 2.5 dimensional hydrodynamical formulation of a disc-outflow cou- 
pling system with a set of self-similar solutions fGhosh & Mukhopadhyay 2009 1, 
we have studied the energetics of the system. The self-similar approach is not 
just a mere tool to solve the non-trivial coupled partial differential equations. The 
accretion flow indeed may exhibit self-similar behaviour. It was shown earlier 
that the black hole system GRS 1915-1-105 is chaotic in nature [IMisra et al. 200411 
that supports the idea of inner disc instability and then turbulence. It was also 
found that the corresponding correlation/fractal dimension to be similar to that in 
the Lorenz system which is a model example of an ideal chaos. The low correla- 
tion/fractal dimension implies possible self-similarity into the system. In reality 
the accretion flows with strong outflows are likely to be geometrically thick and 
advective, which more possibly occur in ultra-luminous and under-luminous ac- 
creting sources. The standard Keplerian model [Shakura & Sunyaev 1973 ] fails to 
describe these two opposite traits of observational signatures, which are possible 
sources of powerful outflows and jets. 

The energetic profiles strongly depend on a and /. They also reveal that out- 
flows and jets are more probable and powerful for super-Eddington accretion flows 
compared to the sub-critical ones, which correspond to ultra-luminous and under- 
luminous sources respectively. For super-Eddington flows (Fig.l), the disc gets 
truncated at larger radii compared to the sub-Eddington ones. This means that the 
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outflow/jet may occur at the inner region of the disc and the disc-corona transition 
region may form in the vicinity to the central star in case of sub-Eddington flows. 
However, a slight departure of y from ~ 1 .4 to a higher value leads to an abrupt 
decrease in the energetics of the flow. We have shown that at super-Eddington ac- 
cretion flows (M ~ 10^ M a ), with an appropriate choice of a and /, the luminosity 
calculated from our model is L ~ lO'*' erg/s corresponding to a black hole of mass 
M ~ IOMq. Such a high luminosity is commonly observed in exotic ULXs in 
external galaxies or in our own super-critical accretor SS433. 

With the appropriate choice of free parameters a, f, M and M, we have shown 
that our model can reproduce the luminosities observed in ultra luminous sources 
and other family of high luminous AGNs (Fig. Ic), which probably accrete 
super-critically. Similarly, the luminosity observed in under-luminous AGNs (e.g. 
Sagittarius A*) can also be reproduced by our model at / — > 1 (Fig. 3c). This 
particular situation is similar to ADAF ^Narayan & Yi 1994| . The limitation of 
ADAF can be addressed by incorporating outflow from the disc (as of Ghosh & 
Mukhopadhyay 2009), which may turn the disc into a non-radiative accretion flow 
and thus can explain luminosity profile of under-luminous sources. As our solu- 
tions can explain the observed luminosities of both ultra-luminous sources and 
under-luminous AGNs with the appropriate choice of free parameters, we should 
put our model to more observational tests. Although we have neglected the contri- 
bution of the magnetic field, whose importance for super-critical accretion flows 
is probably not significant [ |Ghosh & Mukhopadhyay 2009[ , the inclusion of the 
magnetic field will render the system to be more realistic. Such a work will be 
pursued in the course of time. In the follow up work, we plan to analyse a full 
scale numerical solution to conform/verify the results discussed here. 
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Figure 1: Variation of (a) mass outflow rate, (b) outflow/jet power, (c) luminosity 
in unit of Mc^, (d) Bernoulli's constant, as functions of radial coordinate for super- 
Eddington accretion flows, (e) Variatiop^of Bernoulli's constant as a function of 
z. Solid, dotted, dashed curves are for a = 0.3 and / = 0.4, 0.5, 0.7 respectively. 
Long-dashed, dot-dashed curves represent flow with a = 0.01 and / = 0.4,0.5 
respectively. Other parameters arey = 1.4, corresponding ;S ~ 1/3, andz = 5. 
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Figure 2: A comparison of various quantities between two different y for super- 
Eddington accretion flows as functions of radial coordinate. Solid and dotted 
curves are for y = 1.4 and y = 1.44 respectively. Other parameters are or = 
0.01,/ = 0.4,z = 5. 
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Figure 3: Same as Fig. 1, but for sub-Eddington accretion flows. Solid and dotted 
curves are for a = 0.3, / = 1 and a = 0.01, / = 0.9 respectively. Other parameters 
are 7 = 1.6, corresponding jS ~ 0.89, z 



